Microscopic theory for interface fluctuations in binary liquid mixtures 
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Abstract 

Thermally excited capillary waves at fluid interfaces in binary liquid mixtures exhibit simultaneously both density and 
composition fluctuations. Based on a density functional theory for inhomogeneous binary liquid mixtures we derive an effective 
wavelength dependent Hamiltonian for fluid interfaces in these systems beyond the standard capillary-wave model. Explicit 
expressions are obtained for the surface tension, the bending rigidities, and the coupling constants of compositional capillary 
waves in terms of the profiles of the two number densities characterizing the mixture. These results lead to predictions for 
grazing-incidence x-ray scattering experiments at such interfaces. 
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I. INTRODUCTION 



If two thermodynamically coexisting fluid phases are 
brought into spatial contact via suitable boundary con- 
ditions, an interface forms which interpolates smoothly 
between the bulk properties of the coexisting phases. 
For more than a hundred years substantial theoretical 
and experimental efforts have been devoted to resolve the 
structural properties of this transition region (see, e.g., 
Refs. 0,0). 

The reason for the persistence of these challenge resides 
in the difficulty to describe the simultaneous occurrence 
of bulk fluctuations reaching the interface and of capil- 
lary wave-like fluctuations of the local interface position 
For the simplest example, i.e., the liquid- vapor inter- 
face of a one-component fluid, the concept of an effective 
interface Hamiltonian leads to quantitativepredictions 
for a wavelength-dependent surface tension |3, |j, la, la, III > 
which has been confirmed experimentally for the interface 
structure factor down to microscopic length scales for a 
wide variety of one-component fluids H, M EH d E3] • 
For these systems the wavelength dependent surface ten- 
sion j(q = 27r/A) is a function of temperature T and a 
functional of the interaction pair potential of the fluids 
particles. The macroscopic surface tension 7 = 7(0) of 
the liquid-vapor interface is obtained for q — > whereas 
j(q) decreases for increasing values of q, reaches a sub- 
stantial minimum, and increases again for large q. This 
decrease of 7(g) is in a accordance with simulation data 
13 0, E3 which, however, have not yet confirmed the 
predicted and experimentally observed re-increase of "f(q) 
at large q. 

The present work aims at extending this analysis to the 
case of binary liquid mixtures composed of species A and 



B. This is motivated by the following reasons: 

(i) Binary liquid mixtures are governed by three pair 
potentials «%(r) for the A- A and B-B interac- 
tion between the like species and the A-B inter- 
action between unlike species. Provided a wave- 
length dependent surface tension can be introduced 
analogous to the one for one-component fluids, it 
will therefore be a functional of three pair poten- 
tials. By exchanging systematically one of the 
two components by a sequence of molecules with 
a quasi-continuously changing architecture, this 
might open the possibility to tune the shape of 
the function 7(g) and thus to create new interfa- 
cial phenomena. 

(ii) Whereas for one-component fluids two-phase coex- 
istence is confined to a liquid- vapor coexistence line 
described by the chemical potential fi (T), in bi- 
nary liquid mixtures two fluid phases can coexist 
on a two-dimensional sheet in their thermodynamic 
parameter space (ha ,Hb,T) spanned by the chem- 
ical potentials HA an d Hb of the two species and 
temperature (see Fig.0. This allows one to vary 
the thermodynamic state of the system over a con- 
siderably larger parameter space without loosing 
two-phase coexistence, which in turn increases the 
possibilities to vary 7(g) by changing thermody- 
namic variables such as the composition. 

(iii) Generically, for one-component systems liquid and 
vapor are the only fluid phases and thus liquid- 
vapor interfaces are the only possible fluid inter- 
faces in such systems. Binary liquid mixtures ex- 
hibit various fluid phases: a mixed vapor phase, a 
mixed liquid phase, an A-rich liquid phase, and a 
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B-rich liquid phase, separated from each other by 
sheets of first-order phase transitions which inter- 
sect along a triple line of three-phase coexistence 
and which are delimited each by lines of critical 
points (see, e.g., Refs. [H [H [Hf and Fig.Q}. De- 
pending on the relative integrated strength of the 
attractive parts of the aforementioned three pair 
potentials, there is a wide range of rather differ- 
ent topologies of the bulk phase diagrams of binary 
liquid mixtures 0]. These topologies of the bulk 
phase diagrams allow for four distinct types of fluid 
interfaces: vapor|mixed fluid, vapor|A-rich fluid, 
vapor | B- rich fluid, A-rich liquid|£>-rich liquid. In 
contrast to one-component systems this offers the 
possibility to vary significant features of fluid inter- 
faces without changing the underlying interaction 
potentials but only the thermodynamic state. 

(iv) The description of inhomogeneous binary liquid 
mixtures requires two number density profiles, 
Pa{z) and ps(z), where z denotes the distance 
from the mean interface position along the z axis. 
In many cases it is suitable to introduce instead 
the total number density p(z) = pa{z) + Pb{z) 
and the concentration X(z) = pa{z) — Pb{z) as 
linear combinations. Whereas for one-component 
systems it is straightforward (as in Ref. |3|) to as- 
sign a local liquid-vapor interface position f(x,y) 
to a given density configuration as the position of 
an isodensity surface (e.g., points s where p(s) = 
(puq + Ptiap)/2), such a construction is not clear 
from the outset in the presence of two fluctuating 
densities. Thus, the study of fluid interfaces in bi- 
nary liquid mixtures raises the challenging concep- 
tual issue how and to which extent they can be 
described microscopically in terms of an effective 
Hamiltonian TL[f] and a wavelength dependent sur- 
face tension j(q). 

(v) It requires special care to prepare a bona fide 
one-component fluid. Naturally, systems come 
as multicomponent samples. Generally, segrega- 
tion phenomena occur at their interfaces, which 
might influence significantly the interface fluctua- 
tions and vice versa. By choosing suitable series of 
molecules of related architecture and appropriate 
concentrations, binary fluid mixtures offer the pos- 
sibility to interpolate systematically between the 
material properties of the corresponding limiting 
one-component systems, which generates substan- 
tial application perspectives. Finally, binary liq- 
uid mixtures can serve as rudimentary polydisperse 
systems as they occur in colloid suspensions. The 
study of interfacial properties in such systems has 



become very rewarding because they can be ana- 
lyzed in great detail by direct optical techniques 
|2(1 l2l| , allowing for quantitative comparisons with 
theoretical predictions on the scale of the particles. 

As mentioned above, two types of fluctuations occur si- 
multaneously at interfaces: (a) fluctuations of the den- 
sity as they occur in the bulk on length scales up to the 
bulk correlation length £ ; (b) in the absence of gravity 
and for large system sizes the mean position of the in- 
terface can be shifted without cost of free energy. This 
gives rise to thermally excited Goldstone modes lead- 
ing to lateral fluctuations of the local interface position, 
with wavelengths reaching macroscopic values. Depend- 
ing on which type of fluctuation is emphasized, originally 
two different approaches for the theoretical understand- 
ing have emerged. 

As put forward by van der Waals 22], the first ap- 
proach leads to a so-called intrinsic density profile which 
interpolates smoothly between the constant densities in 
the coexisting bulk phases. The interface is laterally flat 
and is kept in place by boundary conditions or a small 
gravitational field acting alongthe interface normal. The 
width of the intrinsic profile 0, HjJ 0] is given by the 
bulk correlation length, which diverges upon approach- 
ing the critical point T c of the corresponding two-phase 
coexistence, reflecting the disappearance of the interface 
at T c . Accordingly, the van der Waals picture is expected 
to capture the interfacial properties at elevated temper- 
atures close to T c . 

The second approach, conceived by Buff, Lovett, and 
Stillinger ji^, describes the width of an interface as a 
result of capillary-wave like fluctuations of a step-like in- 
trinsic density profile. Here only the local interface po- 
sitions are the statistical variables. The resulting mean 
density profile attains the bulk values like a Gaussian 
whereas the van der Waals approach yields an exponen- 
tial decay for short-ranged forces between the fluid par- 
ticles or inverse power laws in the presence of dispersion 
forces 13 . Within the capillary- wave model the width of 
the mean interface diverges upon switching off gravity or 
increasing the lateral system size. This roughening effect 
is missed by all available van der Waals approaches. On 
the other hand the capillary wave model misses the fact 
that the interfacial width diverges upon approaching T c 
on the scale of the bulk correlation length £. 

Accordingly one can state that the van der Waals ap- 
proach captures fluctuations on the length scale of £ and 
below and is suitable at high temperatures whereas the 
capillary wave approach is valid at low temperatures and 
captures the fluctuations with wavelengths larger than 
£. In Ref. Q these two approaches have been reconciled 
by considering intrinsic density profiles, as obtained from 
density functional theory, which undergo fluctuations of 
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their lateral positions. Density functional theory provides 
expressions for the cost in free energy of such density con- 
figurations relative to the free energy of a flat interface. 
This yields an effective interface Hamiltonian H[f] and 
thus provides the statistical weight of interfacial fluctua- 
tions /. This statistical weight can also be used to cal- 
culate correlations of the local interface normals [13, 13 . 

Inspired by the motivation described above, the 
present work extends the concept of Ref. |3| to the de- 
scription of fluid interfaces of binary liquid mixtures. Af- 
ter a brief discussion of the bulk phase diagrams of binary 
liquid mixtures (Subsec. HI A|) we introduce the density 
functional theory which we use as the starting point for 
the description of spatially inhomogcneous fluids (Sub- 
sec. ^TOJ. We define the effective interface Hamiltonian 
H for mixtures in Subsec. Hl CI After discussing the crude 
approximation of steplike density profiles fSubsec. lIID|l . 
in Subsec fll El we introduce the central approximation 
which we actually use for further calculations. It in- 
volves the influence of the curvatures of the iso-density 
contours on the density profiles which has turned out 
to be crucial in order to describe the fluctuations of a 
liquid- vapor interface (see Ref. Q). Since this approach 
cannot simply be transferred to binary liquid mixtures 
requiring two iso-density contours, Sec. ^ is closed by re- 
marks about how to overcome these additional problems. 
In Subsec. lll Fl and Appendices^] and D~D we present the 
explicit expressions for the effective interface Hamilto- 
nian Ti. based on the above-mentioned approximations. 
In order to be able to make predictions for scattering ex- 
periments from such interfaces, in Sec. lIIII we present a 
Gaussian approximation of the effective interface Hamil- 
tonian Ti using various representations and we discuss 
the resulting contributions (Subsecs. lIII AI O. In Sec.HVl 
we analyze the temperature and the composition depen- 
dence of structural properties of interfaces in binary liq- 
uid mixtures as inferred from the correlation functions. 
We summarize our results in Scc.lVl 



II. EFFECTIVE INTERFACE HAMILTONIAN 

In this section we derive an effective interface Hamil- 
tonian Ti. for the interface between two fluid phases of 
a simple binary liquid mixture consisting of spherical 
particles with radially symmetric interaction potentials. 
The system with its interface is described microscopi- 
cally in terms of a simple, but for the present purpose 
appropriate grand canonical density functional. For 
each of the two equilibrium particle density distributions 
we specify implicitly an iso-density contour as its inter- 
face surface assuming that this captures the interface 
structure of the mixture as a whole. The interface 



V-A+V-B 



vapor 



B rich 
liquid 




rich 
liquid 



Fig. 1: Schematic bulk phase diagram of a binary liquid mix- 
ture in the thermodynamic parameter space (/ia + /is , HA — 
fiB,T) spanned by the chemical potentials \xa and \xb of the 
two species and temperature T. As mentioned in the intro- 
duction, the phase diagram exhibits a vapor phase, a mixed 
liquid phase, an A-rich liquid phase, and a B-rich liquid phase 
separated by sheets S\ and S2 of first order-phase transitions 
which are bounded by lines L\ and L2, respectively, of second- 
order phase transitions. The line of intersection between Si 
and S2 represents the triple line TL where three phases coex- 
ist. T ccp denotes the critical end point. In the present context 
of fluid interfaces the solid phases of a binary liquid mixture, 
which occur at high pressures (~ {ia + Mb) or l° w temper- 
atures, are omitted for reasons of simplicity. Further details 
can be found in Refs. |16L llTL llSIl . 



effective Hamiltonian, which counts the cost in free 
energy to deform the interfaces from a given reference 
configuration, is defined as the difference between 
two grand canonical potentials corresponding to two 
different surface configurations. Further simplifications 
are made to express this Hamiltonian explicitly, rather 
in terms of the surfaces, in terms of the yet unknown, 
inhomogeneous densities. Thus, by construction, the 
microscopic interactions between the particles are taken 
into account transparently, which finally lead to effective 
interactions between the surfaces. To a large extent the 
functional dependence on the interaction potentials is 
kept general. Ultimately, for numerical evaluations, we 
assume long-ranged attractive dispersion forces. 

The normal of the mean interface of the binary 
liquid mixture is taken to be oriented along the z-axis 
such that, for instance, the liquid phase and the vapor 
phase of the mixture are approached for z — > —L and 
z — > +L, respectively (see Fig. 0). Density functional 
theory assigns a free energy to each density configuration 
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Fig. 2: Sketch of the liquid-vapor interface region of a binary 
liquid mixture with its normal oriented along the z-axis; for 
each of the two fluctuating densities of the mixture an inter- 
face position can be defined (dark and light grey, see the main 
text, and Eqs. @ and JSJ). In order to derive the interface 
Hamiltonian, two interface configurations are considered: the 
reference surface configuration is given by constant iso-density 
contours for the two species A and B at z — c\ and z = C2, 
respectively, while a non-flat configuration with iso-density 
contours z — fi(x,y) and z = fiix,])) varies laterally around 
the position ci and C2, respectively, so that (/i) = ci and 
(/b) = C2. Additionally, the non-flat surfaces are assumed 
to fluctuate mildly without overhangs or bubbles, so that a 
Monge representation with fi(x,y), i £ {1,2}, can be used. 
Vapor and liquid should be understood as representations of 
two fluid phases including a mixed liquid, an yl-rich liquid, 
and a B-rich liquid (see Subsec. Hl Al . 

such that the equilibrium configuration minimizes the 
functional and yields the corresponding grand canonical 
potential. As the natural reference configuration we 
choose what we call the flat state, in which the iso- 
density contours are laterally constant surfaces and do 
not vary with R = (x,y) (see Fig.[2}. If present, gravity 
points into the negative z-direction. 



A. Bulk phase diagram of binary liquid mixtures 

As stated in the introduction, binary liquid mixtures 
are composed of two species, called A and B particles. At 
high temperatures these particles mix in a gaseous phase. 
Upon lowering the temperature the mixture exhibits a 
phase separation into a gas phase of low density and a 
liquid phase of high density. In Fig.^ this phase separa- 
tion is indicated by the sheet Si with pa + Pb as a mea- 
sure of the total pressure of the system. At sufficiently 
high temperatures in both these phases the two species 
remain mixed. A further decrease of the temperature 
leads to an additional phase separation of the fluid phase 
into an A-rich liquid phase and a B-rich liquid phase (see 
sheet 5*2 in Fig.[Tjl . In the following any pair of the mixed 
gas, mixed fluid, A-rich liquid, and £?-rich liquid are de- 
noted as liquid and vapor. Their coexistence corresponds 



to a point on S\ or S2 and, for instance, an increase of 
temperature at coexistence delineates a path on Si or S2 
approaching the line of critical points L\ or L2, respec- 
tively. On the other hand, changing the composition of 
the mixture at a fixed temperature at coexistence corre- 
sponds to a path on S\ or S2 intersecting a horizontal 
(/j,a + Pb,Pa — /^s)-plane in Fig.^ In Sec. lIVI we shall 
discuss our results in two respects: first, we study the 
influence of temperature and, second, we shall keep the 
temperature fixed and consider composition variations. 

B. Density Functional Theory 

We consider a grand canonical density functional for a 
two-component fluid which consists of particles A and 
B with a spherically symmetric interaction potentials 
Wij(l r l) = Wji(l r l)i where the indices i,j £ {1, 2} refer 
to the species A and B. Following standard procedure 
|29| the interaction potential is split into a short-ranged 
repulsive part Sij(r) and an attractive long-ranged part 
Wij (r) . For a system of volume V = 2L A, where A is the 
(flat) interfacial area and 2L is the macroscopically large 
extension in z direction, a simple version of the grand 
canonical density functional Q reads: 

n[ Pl (r),p 2 (r)] = ^ s [ Pl (r) )P2 (r)] 

+ ^ [d 3 r (w + vrWMr) 
<=i Jv 

2 

+ \tl I d " r l dV «%(|r-r|)Mr) Pj (r'). (1) 
2 iJ= i Jv Jv 

Here, Pi(r) is the number density of the particles of 
species i £ {1,2} at r = (x,y,z) = (R, 2), and 
J- bs [pi(r), P2(r)] is the reference free energy functional 
of a system governed by the short-ranged contribution 
Si } ;(|r|), expressed suitably in terms of a hard-sphere sys- 
tem. In the following, we use J- bs [pi(r), /92(f)] within a 
local density approximation: 

^ h5 [ Pl (r), P 2(r)] - f d 3 r h( Pl (r), p 2 (r)) . (2) 

In Eq. the chemical potential of species i is denoted by 
/ij, while V^lr) represents its external potential, which 
in our case will be gravity acting along the negative z- 
axis. The attractive part of the pair interactions is given 
by tUy(|r|) = Wij(r). To a large extend our reasoning 
will not depend on specific choices for h(pi(r), p%(r)), 
Vf x , and Wij(r). This will be required only for quanti- 
tative presentations. Actually Wij(r) should be replaced 

(2) 

by the direct correlation function c£. (r) which, however, 
reduces to Wij (r) for large r. This replacement also does 
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not alter our main results. 

With the notation R = (x,y) we introduce the bulk den- 
sities 



P 



(3) 



characterizing the vapor (pf) and the liquid phase (p~) 
in the general sense described above. In order to de- 
scribe density configurations as shown in Fig. we intro- 
duce p Ci {z) and Pfi(j) as the density profiles of species i 
which take a fixed value pi at the position z = a for a flat 
configuration and at z = /s(R) for a non-flat configura- 
tions, respectively. For the non-flat iso-density surfaces 
we assume a Monge parameterization (see Fig-EJl ■ Thus, 
the crossing criterions are 



and 



Pa(z = Ci) = p t 

Pf^ z = M R )) = Pi 



(4) 
(5) 



The indices c, and /, indicate that these functions of z 
only and of r = (R, z) take the constant value pi at z = Ci 
and at Z = /i(R), respectively. Reasonable choices for pi 
would be pi := {pj + pf)/2 or the analogue of the Gibbs 
dividing surface concept in the one-component fluids (see 
also Fig.|3|); however, our results do not depend explicitly 
on the choices of pi. Finally we introduce the density 
differences 



Pi 



(6) 



In the following we choose p Ci (z) and p/ ( (r) such that 
they minimize Eq. Q under the constraint given by 
Eq. J3J and Eq. J5J, respectively: 



Sn[ Pl (r),p 2 (r)} 
5pi{r) 



0, i£ {1,2}, 



(7) 



c,f 



where 



denotes the functional 



(dCl[pi,p 2 }/5pi(r) 

derivative of Q w.r.t the density pi, under the constraint 
c (see Eq. (@J) or / (see Eq. (jSJ), respectively. 
Within density functional theory, p Ci {z) and p/ ; (r) are 
equilibrium density profiles in the sense described before. 
Inter alia, Eq. will allow us to eliminate the explicit 
dependences on the chemical potentials pi in our analytic 
expressions; for this purpose it is sufficient to use Eq. (jTJ 
only for the profiles p Ci (z). Without constraint Eqs. 
and J2J) lead to 



5n[p 1 {v) lP2 {v)] 



9 pi /i( Pl (r),p 2 (r))+ Mi + y-(r) 

2 

/ d 3 s w rt (\v-s\) Pj (s). (8) 

3 = 1 JV 



Up to here there is no construction scheme provided for 
determining p Cl and p C2 . One can take solutions of Eq. (JHJ 
for V? x = and shift the pair such that, e.g., the condi- 
tion p Cl (z = c\) = pi is fulfilled (see Eq. (@}), but in gen- 
eral p C2 {z) will not have the property p C2 (z — c 2 ) = p~2- 
This shows that there is only one degree of freedom in 
shifting, i.e., c\ — c 2 . Thus, c 2 is not a free parameter but 
depends on c\ , which means that f2 cannot be minimized 
for arbitrary pairs (ci, c 2 ). As a consequence the effective 
interface Hamiltonian TL depends only on the difference 
ci — c 2 , but we shall treat c 2 formally as a free param- 
eter, which indicates the position of the planar interface 
ofp C2 (z). 

The free energy density h(^pi,p 2 ^j = h' d (pi, p 2 ^j + 
h ex (pi, p 2 ) of the hard sphere part consists of an ideal gas 
contribution ft, ld (pi,/52) and an excess part h ex (p\, p^. 
Our analytical formulae derived below will not depend on 
its functional form; for numerical calculations, however, 
the Carnahan-Starling ex pres sion /i cs (pi,p2) for the ex- 
cess contribution is used |30|. With (/3 = (ksT)^ 1 and 
the thermal de Broglie wavelength A t h, this means 

2 

h' d (p 1:P2 ) - /r 1 J> (ln(A t 3 MA ) - 1) (9) 



and 



h CS {pi,p 2 ) = - n a In (1 - n 3 ) + - 



m n 2 



n3 



n 3 + (1 - n 3 f In (l - n 3 ) 
367™^! — n 3 ) 2 



(10) 



where 



2 

n a = y^Pfc 

k=l 



Til 



fe=i 



(ii) 



47T 



"2 



The weighted densities n a , a £ {0, 1, 2, 3}, are composed 

(i) 

of the densities pi and the particle radii r of species 

ie {1,2}. 

In order to model the van der Waals forces of simple fluids 
we take for the attractive part of the interactions 

which gives the correct large distance behavior Wij (r 3> 
ro^) — —Wo'r~ 6 . The quantity Wo represents the 

(ii) (i) (i) 

depth of the potential, while r = r Q + r is defined 
as the sum of the particle radii. The functional form of 
Wij (r) for small r is chosen for analytic convenience; most 



5 




9 c (z) 



Fig. 3: Sketch of the planar number density profile p c (z); the 
width 2£ = £ + + £~ of the transition region is roughly the sum 
of the bulk correlation lengths £ + and £~ of the two coexisting 
phases which, in general, differ from each other (see also the 
remarks related to, c.f., Eq. H56[l 1. The length £ is also called 
the interfacial width of the (planar) surface located at z — c 
(see Eq. (0J) . The marked regions are frequently used to define 
the interface position via Gibbs' zero-adsorption criterion: the 
number of particles within the dark and the light domain has 
to be the same for the interface at z = c 1,] . Within the sharp 
kink approximation one has £ = so that a step-like profile 
results (see Eq. 1261 1. 



of our results do not depend on this choice. 

Independent of the explicit form of the potentials twy 
we introduce their integrals 



w$(\B.\,z) 



Scij 



dz i%(|R|, \z\) 



and 



w 



(2) 



|R|,z) := f dzw^(\R\,\z\) 

J Sen 



(13) 



(14) 



where Sc^ := Ci — Cj. They fulfill the symmetry relations 



^ ) (|RU) = (-1)^(|R|,-^), feG{l,2}. (15) 

In the present context the binary liquid mixture is ex- 
posed to a gravitational field acting along the z direction: 



V t ex (z) = Gm z (z- Ci ), 



(16) 



where G is the acceleration of gravity, is the parti- 
cle mass and c, the equilibrium flat interface position of 
species i (see Eq. (0J). In the following the first integral 
of Vf x (z) is frequently used: 



V t {1 \z) := \Gmdz-ctf 



(17) 



In the following subsection the effective interface 
Hamiltonian will be defined on the basis of the density 
functional £1 introduced above. 



C. Effective Interface Hamiltonian 

The envisaged effective interface Hamiltonian TL for a 
binary liquid mixture provides the cost in free energy to 
maintain interface configurations described by the iso- 
density contours /i(R), i € 1,2, relative to flat configu- 
rations Pa(z)- We expect that TL depends on the differ- 
ences /f (R) := /i(R) — Ci. Therefore, we introduce the 
abbreviations 



f(R,R) 



A C (R) 
/l(R') 



(18) 



f(R) := f(R,R), <5/£(R,R') := /f(R) - /?(R'), 
%(R,R') := /i(R) - /i(R')j and 5/ y (R) := /,(R) - 
/j(R). In terms of these quantities, the effective inter- 
face Hamiltonian TL is defined as the difference of the 
corresponding grand canonical potentials: 

H[f (R, R')] := 0[PA (r), P h (r')] - n[p 0l (z),p C2 (z')} . 

(19) 

Our main goal is to derive an explicit expression of TL in 
terms of /f (R). 

We rewrite TL by carrying out partial integrations such 
that TL is expressed mostly in terms of derivatives of 
profiles which are mainly confined to the interfacial re- 
gion and vanish for z — > ±L. Due to Eq. © one has 
|/f (R)| <C L for all R, i.e., the interface deviations are 
much smaller than the sample size. According to the 
structure of f2 in Eq. (JJJ, TL is the sum of four terms. 
The first term, which we shall treat later in Subsec. lll Fl 
is given directly as 



H h (f(R)) 



dz 



-L 



h {Ph ( r )> Ph ( r ))~ h {Pci 0), Pc 2 0)) 

(20) 

The second expression stems from the external potential 
V? x (z) and has the form 

H v (f(R)) := fdzV^iz) [d zPfi (r) - d z p Ci (z)] . 

. (21) 

The third contribution involves the chemical potentials /ij 
and additional boundary contributions which arise from 
the interaction potentials. With the constants 



Pk 



+ 4tt^ [d 2 R[p j w$(\R\,z = L) 

-pJw^(\B.\,z = 0)] (22) 
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it reads 



i=l 



Hfc(f(R)) :=-£ Ki / dz(z- Cl ) [5 2/9/i (r)-5 zPci (z)] 



(23) 

Finally, the contribution to Ti due to the attractive part 
of the interactions can be expressed as 

H w (f(R, R')) := -~ Y] f dz f dz' w ( i 2) (\5R\,5z)x 

A ilj=1 J-L J-L 

'd z p u {v)d z ,p f] {v')-d z p c% (z)d z ,p C] {z')\ (24) 

where <5R := R R and := z - z' . Thus, Eqs. O- 
J2U lead to 



W[f(R)] 



W h (f(R)) +W 6 (f(R)) +W v (f(R)) 



jfjf d 2 Rd 2 R' H w (f (R, R')) 



(25) 



In the following two subsections we analyze two different 
models for the profiles p Ci {z) and p/;(r) in order to ob- 
tain analytic results for Ti. The first approach assumes 
that at the interface position the densities vary discon- 
tinuously between the corresponding bulk values. This 
so-called sharp kink approximation will be discussed in 
Subsec. Hl Dl The second approach fSubsec. Ill E|) is based 
on continuous density profiles and takes the influence of 
the curvature of the iso-density contours on the densities 
into account. 

The validity of our approach is also based on the assump- 
tion that in the thermodynamic limit, i.e., A — > oo, all 
lateral boundary contributions to Ti vanish in Eq. 125(1 . 



D. Sharp Kink Approximation 

The sharp kink approximation replaces the actual 
smooth variations of the density profiles (see Fig.[3J) on 
the scale of the bulk correlation length by step functions: 

Pfi (r) := -ApiBiz-MR^+pr (26) 

where 0(x) is the Heaviside function so that 

d zPh (r) = -A/H6{z-fi{K)). (27) 

Similar expressions hold for p Ci {z) with /i(R) replaced 
by Ci. For one-component fluids this approximation has 
turned out to be surprisingly successful in describing 
liquid- vapor interfaces and wetting phenomena |2(| . 
From Eqs. (|25|l and Eq. I|27|) together with the expansion 
(see Eqs. (|14 ) . (p ) . and flBjl ) 



w^{\SB\,6f ij ^B!))^w iS {\SR\ i Sc ii ) (R, R')] 3 

(28) 



we find 

tt sk [f(R,R')] 

2 

d 2 R ^2<A Pi 



^rn t [/,f(R)f 



MR) x 



Ph 



2=/<(R 

2 



9 Pi h(p cl ,p C2 ) 



d 2 Rd l R 



<>J=1 



Apt Apj 



xm y (|<5R|,<5 Cu )[^.(R,R)] 2 . (29) 

In Eq. 1|29[) the expressions d Pj h\ z= f i — d Pj h\ z=Ci are ba- 
sically not determinable because at least one density has 
to be evaluated at the interface position of the second 
which is unknown. For instance, in order to evaluate 
Pf 2 (R, /i(R)) in the case c\ = c 2 = one would need the 
information whether /i(R) > /z(R) or /i(R) < /2(R): 
the first case yields p/ 2 (R, /i(R)) = pt > the second gives 
P/ 2 (R, /i(R)) = P2- As a consequence, the expressions 
depend on the differences f 2 (R) — f\ (R) and c\ — c 2 and 
even vanish for /2(R) — /i(R) = and c\ — c 2 — 0, be- 
cause each density is evaluated at its isodensity surface 
resulting in the same value (see Eqs. Q, ©). Using the 
expansion 



*=A(R) 



d Pi h(p 1 ,p 2 ) 



J2d Pk d Pz h(p 1 ,p 2 ) 



k=l 



PL 



and similarly for d Pi h(p Cl , p C2/ 



leads to 



Pk 



(30) 



W sk [f(R,R')] 

i=l 



d 2 i? [/p(R)r 



(31) 



h(pi,p 2 ) x 



/ d 2 i? |<5/ y (R)| 2 e( - <y/« (R)) 
JJj 2 Rd 2 R' viijQSRlScy) [<y/P-(R, R) 



Note, that the Heaviside function vanishes if i5cy and 
have the same signs. Its prefactor |/i(R) — /2(R)| 
prevents an appropriate Fourier analysis because the re- 
sulting expressions cannot be ordered in products of fifj, 
where / denotes the Fourier transformed function of / 
(see Eq. (J22J). Therefore, within this sharp kink approx- 
imation, the cost in free energy for deforming the inter- 
face can be studied only for the case /i = f 2 but not for 
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the more general situation /i ^ / 2 . 

For /i(R) = / 2 (R) ee /(R) and c x = c 2 = the afore- 
mentioned problematic expressions in Eq. I|29|l drops out. 
With the Fourier transformation 



vapor 



/(q) := WF" L dnRm)e ^ R - 
/(R) = WF" i d>/(q)e+iqR ' 



(32) 
(33) 



and 



1 

2^ 



d 2 R Wij (\R\,u)e~ lCiR 



/ dRRJ (qR)w i:i (R,z) 
Jo 



(34) 



where J (x) := ^ J e 4a: cost dt is the zeroth order Bessel 
function, 7i sk can be expressed as 



H sR [/(q)] = - / d^|/(q)|" Ge sk + ^7 5k (g) , (35) 

47T J R 2 L J 

with Q sk := Ylj—i A/3j mj and a wavelength-dependent 
surface tension 

i 2 

7 5k (<z) :=-jj; A Pi A Pj [m y (|q|,0)-m y (0,0)] . (36) 

q i,j = l 

Equation (|36l) is the generalization of the corresponding 
result for a one-component fluid 0, 0] assuming a single 
steplike interface in the binary case. For fluids governed 
by dispersion forces (Eq. I |12| l) one obtains in the limit of 
long wavelengths l/q 

7 5k (? - 0) = 1 A <°* A ^ w o j) ( ri o j) ) 4 x (37) 

(l+i (Z 2 (rfe)) 2 [log(^) + c])+0( (? 4 ) 
with Euler's constant C £ = 0.5772 . . . and C = C e - §; 

1 2 

7o 5k := 7 sk (? = 0) = - ]T A Pi A Pj «,<«) (r^)) 4 

(38) 

is the macroscopic surface tension within the sharp kink 
approximation. At short wavelengths, i.e., q — > oo, one 
finds 7 sk (<? — * oo) — > 0, which means that distortions 
with short wavelengths are insufficiently suppressed. 

While the previous calculations are based on intrinsic 
steplike density profiles, in the following subsection we 
consider the more realistic case of smoothly varying in- 
trinsic profiles, including changes of their shape due to 
local curvatures of their interfaces. 




Fig. 4: Sketch of an interface /(S), S € R 2 , and its normal 
coordinate system (NCS) consisting of a point s G R 3 on /, its 
normal vector n, and the normal distance u. Thus each point 
P has two representations, either as a vector r or as s + u n 
within the NCS (see Eq. 1391 '). Additional assumptions are 
required to assure the unique equivalence of both (see main 
text). In the present picture the condition |it| < R m in for all 
u is violated because the normal distance of the point P is 
larger than the radius of curvature of the right bump. 



E. Curvature Expansion 

In this subsection we consider continuous density pro- 
files p Ci ( z ) (see Fig.[3J) and p/ ; (r). The thickness of the 
transition region or the width of the interface is of the 
order of the bulk correlation length £. 

In order to take the influence of local curvatures on the 
density profile pf t into account, first we introduce normal 
coordinates for each surface /i(R) followed by a trans- 
formation of the density p f ( (r) to its normal coordinate 
system. Second, the transformed density is expanded in 
powers of local curvatures. (For the following general 
remarks we omit the index i.) 

To this end we consider the points s(S) = (S, /(S)) of 
the Monge parameterized surface /(S), the normal vector 



n(S), and the map T 



(see Fig.0J so that 



T(S,u) := s(S)+wi(S) . 



(39) 



Thus, each spatial point r can be expressed in terms of a 
point on the surface s and its normal distance u from the 
surface. However, finding S and u for a given point r, i.e., 
finding the solution of the equation r = T(S(r), u(r)), is 
generally not a trivial task. However, in order to obtain a 
unique map T _1 we have to restrict the range of values of 
u to {—Rmin, Rmin) where R m in > denotes the absolute 
value of the minimal radius of curvature of the manifold 
/(S). Therefore, the constraint \u\ < R m in guarantees, 
that the Jacobian J/ of the transformation T, 



J/(S,u) = ^g{S)(l-2H(S)u + K{S)u 2 ), (40) 

does not vanish in the domain R 2 x (— R m ini Rmin) ■ 
Here, H = H(S) is the local mean curvature, K = K(S) 
is the local Gaussian curvature, and g = g(S) is the met- 
ric of the manifold /(S) which in Monge representation 
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takes the form g(S) = 1 + (V/(S)) 2 . 

If /(S) is an iso-density contour of p/(r), we may 
write for points (S,w) el 2 x (— R m in, Rmin) 



Pf( r ) = P/(S,u). 



(41) 



These expressions hold also for a fiat surface which 
results in p c {z) — p c {u) with u = z — c. 

Similar to Eq. we assume 

p fi {S,u^±R min ) = pf , i e{l,2}. (42) 

Since Eq. i|42[) can be strictly valid only for macroscopicly 
large values of u, we assume that R m in is sufficiently 
large so that Eq. (|4*2"|) is fulfilled for all practical purposes. 
Now, we propose an expansion of the transformed density 
profile pf i (S, u) into powers of the local curvatures Hi = 
Hi{S) and Ki = Ki(S), i £ {1,2}; 



/5/,(S,u) 
<Wi(S,u) 



Pa{u) + Sp fi {S,u) 



a, 13=0 



{2H l TKfp HaK ,{u) 



(43a) 
(43b) 



For each term this implies a factorization of the depen- 
dencies on the lateral coordinates S and the normal dis- 
tance u, reflecting the condition that the width of the 
interface £ should be small compared with the minimal 
radius of curvature, i.e., £ <C Rmin- For the following 
calculations, it is not necessary to specify the functions 
p\(u), A £ {H,K,H 2 ,HK, . . .}, which depend only on 
the normal distance u but are so far unknown explicitly. 
However, for quantitative predictions one has to use a 
model for p\{u) (see, c.f., Subsec lTlI A(l . 



F. Mean Surface Approximation 

Except for H.h (see Eq. (120(1 ) the formulas derived 
above can be used to transform and to expand the var- 
ious contributions of the Hamiltonian Ti.. For Tih both 
densities have to be evaluated at the same spatial point, 
but there is no rule telling which normal coordinate set 
should be used for the transformation, i.e., which local 
curvatures have to be used. In order to resolve this is- 
sue we construct a mean density distribution p* (r) with 
a corresponding iso-density contour /*(R), such that 
/i(pi(r), /?2(r)) = h*(p*(r)) where h* is a function to be 
determined, and use the normal coordinate system as- 
sociated with f*(R) in order to transform Tih and to 
expand h*(p*(r)). Since we know the relation between 
/*(R) and fi(R) explicitly, we are able to express the 
results in terms of fi. We stress that this problem would 
equally arise for more sophisticated density functionals 



beyond the local density approximation used in Eq. J2J). 
We start our approach with an implicit definition of the 
points r* fulfilling the constraint p*(r*) = const. First, 
we assume that the flat configuration can be written as 
(see Eqs. ©-©) 



Pc, (u) = - 4p p(u/£i) + Pi 



(44) 



with an odd function p(— x) = — p(x) so that p(0) = 0. 
Second, we define, with a not yet specified prefactor Ap*, 
the mean density 



/„„V ^^-"" . (45) 

i=l 



The equation for the corresponding iso-density contour 
reads 



p*(r*) 



p ^ Api 

i=i y% 



(46) 



Using the expansion introduced in Eq. (|43 bfl . in lowest 
order this leads to 



A p Ci ( Ui (r*)) 



(47) 



and hence to the condition (using Eq. I|44|l ) with u* 
Ui{r*) 



Mi Un 

4i ?2 



(48) 



Thus, we postulate that the normal distances u* between 
fi and a point r* on the iso-density manifold /*, mea- 
sured in units of the width £j of the corresponding inter- 
face, are equal (see Fig.[5J); this is a construction scheme 
for /*. In general, the lateral coordinates and S^ 2 ' 
which belong to the same point r* are different. We write 
Si(r*) := (SW(r*),/i(SW(r*))) = s* for the correspond- 
ing point on the surface fi while (r* ) = n* denotes its 
normal vector there. Expressing r* as 



ai s x + d2 s 2 + Q.\ n 1 + a2 n 2 



(49) 



with the coefficients a,i,cti £ R to be determined. In 
combination with Eq. (|48|l one obtains 



= s - I — : — n 



6- 







(50) 



with a x = Ci/(6+£i n^n|), a 2 = Ca/(6+6 nfn|), and 
coefficients Ci £ 

For the special case s* = s 2 one has r* = si = s 2 which 
leads to the relation a\ + a 2 = 1 and oti ~ |s* — s 2 |. For 
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Fig. 5: Sketch of the construction scheme for the mean sur- 
face points r*. Since the points r* represent an iso-density 
manifold /* of the mean density p*(r) they are determined 
in terms of the distances u* between the surfaces ft and r* 
(see Eq. 1481 1 assuming that, normal to each surface fi, the 
density profile has the form given by Eq. 1441 . 



symmetry reasons we set a± = a 2 = \ in order to treat 
the surfaces equally. We now consider the case /2(S) = 
-/i(S) and £i = 6, which implies = S< 2 ) = S and 
r* e z — 0, where e z is the unit vector in the z direction 
(see Fig.0). With Ci = % ( s 2 - s i) n 2 and C2 = % (a* - 
s|) n* this leads to 



■/i(S)V/i(S) 




and u\ = -vOT) /i(S) 



(51) 

so that Pythagoras' theorem, uf = (R* - S) 2 + /i(S), 
is fulfilled. (For different choices of Q this is generally 
not the case.) This means that for /2(S) = — /i(S) the 
manifold /* is the plane z — 0. Using the same choice 
for Q in the case 6 ^ 6j Eq. t4"9"|) yields 



6 H - st] n* 2 



2(6+6 n in *) 
so that the distances u* fulfill 



2(6+6ntn*) 



6 



21 

6 



fl'2 
6 



6 



(53) 



Nonetheless, Eq. I|49|) still is an implicit expression for r* 
which represents approximately points on the iso-density 
manifold /* of p*(r). 

Since f\ and $2 are assumed to not exhibit strong vari- 
ations on short scales, this translates to /* so that r* 
allows for a Monge parametrisation r* = (R*,/*(R*)), 
too. Hence in Eq. I|49ll we can use a Taylor expansion 
/i(SW) re /i(R*) + V/i(R*) (S« - R*) which leads in 
lowest order to (Eq. lf52|0 



/*(R*) 



if 

2 ^ 

i=l 



/i(R*) 



, e ■■= 



266 
6 + 6 



(54) 



A more sophisticated calculation, which takes additional 
curvature corrections in Eq. I|47|l into account, i.e., using 
the next higher order terms in Eq. I|43b(l , shows that cor- 
rections to Eq. are of the order 0(f™ /™ ,n + m > 3). 
Since in Sec. lIIII we shall consider the Hamiltonian Ti. 
within a Gaussian approximation the expression in 
Eq. H54|) is sufficient. Furthermore, from Eq. (|54|l it 
follows that the surface with a smaller interfacial width 
£ contributes stronger to /* . 

To summarize Sec.[n] from the density functional 
fl in Eq. Q for a binary liquid mixture, in Eq. 1191) 
we have introduced an effective interface Hamiltonian 
Ti. by specifying an iso-density contour fi for each 
density profile pi as its corresponding interface, which 
compose the interface of the mixture as a whole, and by 
comparing them with the corresponding flat reference 
configurations. In order to express Ti. in terms of the 
manifolds fi we used an expansion of the densities in 
powers of curvatures of /, (see Eq. (|43b|) '). Since the hard 
sphere contribution Tih cannot be treated within this 
approximation, we have constructed an effective mean 
surface /* (see Eq. Q54^). which itself is an iso-density 
contour of a composed density p* , so that the curvature 
expansion can be performed regarding /* and p* . The 
results of the curvature expansion of 7i up to second 
order are presented in the next chapter. Since all expres- 
sions would become rather clumsy without using short 
notations, we shall introduce additional abbreviations in 
order to obtain a clear presentation of the structure of 
the formulas. 



III. GAUSSIAN APPROXIMATION 

In the previous chapter we have illustrated the ba- 
sic ideas and have derived the general expressions which 
arise upon introducing the effective interface Hamilto- 
nian. In this section we carry out the curvature expan- 
sion in Eq. (|43bJ) up to second order. Higher order terms 
are given explicitly in AppendixlX! Since in the following 
the profiles p Ci (z) do no longer occur we drop the tilde 
in p c {u) and write p c (u) instead (see Eqs. lf4"T|l - (|43bJl ). 
First, we provide some numerical aspects which enter 
into the graphical presentation given below. Within the 
Gaussian approximation Ti is determined by the profiles 
p c and ph and the interaction potentials Wij given in 
Eq. I|12|) • For p Ci we use Eq. (14411 with an intrinsic pro- 
file p(x) = tanh(x/2) and, guided by Ref. for we 
choose 



Ph, (u) 



C N r^ 



x p(2 x) 
2ir cosh(a;) ' 



i6' (M) 
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with a dimensionless positive number Cm- Comparing 
this expression with the analogous one, pn , for the one- 
component fluid introduced in Eq. (3.27) in Ref. Q with 
a prefactor Ch, one obtains Cm r — Ch Thus, dif- 
ferent from Ref. 3] here we assume, that the prefactor 

(i) 

Cm To does not vary with temperature. This choice 
here translates into that in Ref. |3J if there one takes 
C H ~ C 1 -> for T -> T c . Therefore, p H J&Pi re- 
mains bounded for all temperatures, so that the curva- 
ture influence characterized by pn { vanishes ~ Api for 
T -> T c . In Ref. 0, p H /&Pi diverges ~ £ for T — > T c , 
so that in that temperature range the influence of the 
curvature may even dominate. Therefore we prefer the 
choice given in Eq. I|55l) A more detailed discussion of a 
possible temperature dependence of Ch can be found in 
Ref. [s3. For reasons of simplicity, in the following we 
consider only the case Cat = 1.0. We emphasize that the 
structural properties found in Ref. Q do not change if 
Ch ~ S,^ 1 instead of being constant. 
While the ratio ri /r 2 ^ of the radii of the particles is 
a free parameter, the temperature dependent correlation 
lengths £j are determined by the bulk correlation func- 
tions. In terms of the total bulk density p the concentra- 
tions pi = Xi p fulfill Xi + X 2 = 1 • From the Ornstein- 
Zcrnike theory for mixtures one has 



i 2 



P 



XT 



V f d\r 2 WlJ {r) 



(56) 



with the isothermal compressibility xt = "(f 2 ) 
which can be expressed as 



T,V> 



1 



XT 



d 2 h{ P X 1 ,pX 2 ) 

dp 2 



2 

-E 



d 3 rX l X ] w ij {r) 



(57) 

For the two coexisting phases liquid and vapor the two 
corresponding total densities p ± lead to different values 
Xj- and thus & = £i(p^)- In the subsequent numerical 
calculations we use & = +^ i ~)/2 (see also the caption 
of Fig.®. 



A. General Expression 

As stated at the beginning of this section we consider 
only contributions to 7i up to 2nd order in the deviations 
ff of the local interface height from the flat configuration. 
With f(q) as the Fourier transform of the vector f(R) 
(see Eqs. (JTHJ and (EU), one has 

1 



W[f(q) 



-in 



dfq f T (q)E(g) f(q) 



with 



E(g) :=G(g)+W(g) 



(58) 



(59) 



Here, the matrix G(q) represents the contributions stem- 
ming from gravity (see Eq. (|fj22(l '). the matrix W(q) cap- 
tures the influence of the attractive interaction potentials 
(see Eq. I |B241 ). and the constant matrix K involves hard 
sphere contributions (see Eq. l|B14(l ). The explicit expres- 
sions for G, W, and K are derived in Appendix^ where 
the equilibrium condition for the planar densities p Ci (u) 
(Eq. JSJ) is frequently used to obtain the final form of 
E(g). 

In order to be able to present our results in a compact 
form we introduce the following abbreviations. For an 
integer n > and arbitrary expressions JL = Ai(u, . . .) 
we define the moments (similar as in Ref. |3j) 



S n [A, 



1 

Ap* J- 



du u n Ai(u, . . .) . 



(60) 



With this notation, the matrix elements of G(q) can be 
expressed as (8 a p is the Kronecker symbol) 

Gn = GmiApi [l-2<7 2 {s o [p Hl } 

2 

+ - 4) Hdu Pck ] + q 2 S 2 [p Hk }))] (61) 

fc=l 

with di — + £2) and similarly for G22 by inter- 

changing the indices 1 <-» 2 and (G12 = G21) 

2 

G12 = 2q 2 d 1 d 2 '^Gm i /\p i {s 2 [d u p Ct ] + q 2 S 2 [p Hz ]j ■ 

4=1 

(62) 

The first part of Eq. (|61|l up to S [phA is identical with 
the corresponding expression in Ref. |3j , and is recovered 
by setting Ap2 = and £2 ~^ 00 (which results in d\ = 
and d 2 = 1), which consequently implies /* = f\ from 
Eq. All further parts in Eq. l(6"T]) arise due to the 

presence of a second interface. This is somewhat surpris- 
ing, because the gravity terms of the density functional 
f2 are diagonal in the densities and thus, one expects 
G(g) to be diagonal w.r.t. the surfaces, too. Actually, 
the additional terms in Eq. (|61|l and Gi2(<?) 7^ emerge 
by applying the equilibrium condition in Eq. JHJ) in order 
to get rid of certain hard-sphere contributions from H.^ 
(Eq. I^OI) ). All remaining hard-sphere contributions are 



captured by the matrix K = (Kjj) 



,ie{i,2} 



with 



du d 2 pzPj h(p Cl (u), p C2 (u)) p H , {u)ph 3 (u) . 



(63) 

K has the form expected as the generalization to two 
components of the analogous term k in Ref. 
The matrix W(q) depends on the pair potentials. Hence, 
for k £ {0,1,2} it is convenient to use the short no- 
tation vjrf' [q, 5cij] := w^\q,u' — u" + Scij) for the 
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Fourier transformed interaction potential (see Eq. 1341) ) 
or the integrals of it (see Eas. (fT3|l and {H}), respec- 
tively. Moreover, similar to Eq. (|60|l for an expression 



Aij = Aij(u' , u" , . . .) we define the moments 



du'du" w\j [q, Sci 3 ] Aij (u', u", 



•) 

(64) 



and the differences 

(q, ■ ■ • ) := (q, ■ • • ) - w W (0, • • • ) , (65) 

but on the lhs we suppress the indices c and the square 
brackets around Aij indicating the functional dependence 
on Aij. In addition we use cj(°\. . .) = u>(. . .) due to 
. .) = Wy(. . .) and similarly &2>W(. . .) = . .). 
Then, the entries of the matrix W(g) can be written as 
(d< = d u > , d" = d u " ) 

Wn(g) = -w(0, d,p Cl d»p C2 ) +Soj(q,d>p Cl d»p cl ) 
+ 2q 2 ^Su(q,p Hl (u')d'>p Cl ) - w(0, p Hl (u')d»p C2 ) 



+ ^2(4- S a )u; {1} (0,u'd,p Cz d„p Cj ) 

2 

+ g 4 [2 X)(^-* i i)w {1 >(0 > «'pH 4 («>'Pc*) 
2 

3=1 



(66) 



and (W12 = W 2 i) 
Wi 2 (g) = Lu(q,d>p Cl d»p C2 ) 



(67) 



; \(l~5 tJ )u(q,p Hi ( U i)d„p Cj ) 

+ 2d 1 d 2 w w (oyap Ci d»p C3 ) 

2<i 1 d 2 £«(0,K , Pff i K)3''%) 



w(t?, p Hi (u')p Hj {u")) 



W 22 (q) is obtained by interchanging the labels 1 <-> 2 
in Wn(g). Again, the result for a single interface is in- 
cluded as a limiting case by setting Ap 2 = 0, d\ = 0, and 
d 2 = 1 (or equivalently £ 2 — > 00) in Eqs. lj6U)l and (Iu7|) . 
which gives W12 = W 22 = and /* = /1 from Eq. (fBljl . 
All additional terms are generated by applying the equi- 
librium condition in Eq. © . Although these expressions 
are useful to determine numerically the matrix elements 
Kij(q), the formal structure of E(g) might be more trans- 
parent in the presentation given in Eq. i|C6(l . 



In order to obtain further insight into the nature of TL it 
is useful to transform E into a diagonal matrix. To this 
end, we define 



Xi(q) := Eu(q)+E 12 (q), 
A+(q) := Ai(g)+A 2 (?), 



and 



A"(«) 



detE(g) 



(68) 
(69) 

(70) 



Furthermore, we define a "mean" surface / + , and a "rel- 
ative" surface / _ via 



/ + (q) := 



2 



M?) 



frt Ai(g) + A 2 (g) 



and 



/"(q) := / 2 c (q) - /r(q) , 



//(q) (71) 



(72) 



the fluctuations of which are decoupled within the Gaus- 
sian approximation, in contrast to /f and /f . This leads 
to 



1 



W"[f(q)] = j- / d 2 q ]T A«(g)|r(q)| 2 . (73) 

This resembles some similarity to the decomposition of 
the two-body problem in classical mechanics. It is impor- 
tant to note that A ± (g) are not the eigenvalues of E(q) 
since the coordinate transformation used and defined by 
Eqs. H71fl and (|72|l is not orthonormal. Moreover, this 
transformation makes sense only for Xj 7^ 0. Therefore, 
the limiting case of a single component is better discussed 
in terms of Eqs. (|58|l and (|59|l as mentioned above. 



B. Energy Density A + (g) 



The quantity A. + (q) can be written as 
A + (q) = Gg+(q)+q 2 j+(q) 



(74) 



where Q + {q) = Qi{q) + Q2 with Qi = rrijApj (l 
2q 2 d [p Hj ]) (Eq. ifBlflV and 7+ (3) takes the form 



7 + 0?) ■= Yl ^io) 



(75) 



with 



8Lj(q,d>p Ct d»p C] ) , 
Jijil) = — 2 — + 26u{q,p Hi d»p Cjj 



<Jj(q,p Hi p Hj ) +K ij + 2ujW(0,p Ht d //Pcj ) . (76) 



This formula is the generalization of the corresponding 
result derived for a one-component fluid except that 
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Fig. 6: Normalized wave-vector dependent surface energy 
7 + (g)/7 + (0) as defined by Eq. (17511 for different temperatures 
T = 0.7... 0.99 T c , where T c is the bulk critical point of the 
binary liquid mixture in coexistence with its vapor, computed 
from the grand canonical density functional in Eq. Q. The 

and 



■ (2) / (1) 

size ratio is To /re, 



„( 22 ) /,„( n ) 



M 11J = 1.05, 



1.001 with r X 12) 
„( 12 ) /,„(") 



' /w„ lx> = 0.5 for the depths of the 
interaction potentials (see Eq. 1121 '). 8c\2 is set to 0, but 
7 + (g)/7 + (0) depends only weakly on <5ci2 for |5ci2| < 5ri 12 '. 
Here, Cjv = 1 (see Eq. 1551 1: another choice Cjv > 1 mod- 
ifies mainly the bending constants Ky (see Eq. (IB141 1 and 
thus the curves increase stronger for larger values of q beyond 
the minimum which is also shifted to larger wavelengths 1/q. 
The opposite behavior is observed for Cjv < 1. For low tem- 
peratures the minimum is rather deep and occurs at large q 
values. Upon raising the temperature the minimum becomes 
more and more shallow and moves to smaller values of q, too. 
This is the same qualitative behavior as observed for the one 
component case |3|. 



the term uj( 2 \0, PHid"p Cj ) shows up additionally. j + (q) 
is plays the role of a wavevector dependent surface ten- 
sion for / + , which then behaves similarly as for a single 
interface. Since according to Eq. (|71|l f + {q) is a linear 
combination of the surfaces, /? can be considered as the 
prime or mean surface of the binary fluid. The functional 
form of 7 + (<z) is shown in Fig.|B|for various temperatures. 



C. Energy Density A (q) 

A~(q) is determined by Eq. I|7UI) . It exhibits a more 
complex structure than A + (q). The explicit expression 
for det E(g) is given by Appendix^] The intrinsic behav- 
ior of the surfaces is given by the pair potentials of the 
particles alone, independent of the external field. Thus 
in the absence of gravity, i.e., for G = 0, one obtains the 
undisturbed energy density A~(q) = A~(q, G = 0) of the 
different surface configurations f~. Similar to A + (q) it 
can be decomposed into a wave- vector dependent surface 
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Fig. 7: Normalized wavevector dependent free energy density 
Ao(g)/Ao(0) as given by Eq. (1771 for various temperatures 
T = 0.7 . . . 0.99 T c . The choices for the interaction potentials 



are r jr = 1.001 for the radii and w /w = 1.05, 
w^, 12 " 1 /!()„ = 0.5 for the interaction strengths (see Eq. 1121 1. 
5ci2 is taken to be 0. Distinct from the normalized surface 
tension 7 + (g) shown in Fig. [(J A„ (g) shows a monotonic in- 
crease at low temperatures, but forms a minimum at a nonzero 
value of q if the temperature is increased. As for 7 + (?), this 
minimum is shifted to larger wavelengths upon further in- 
creasing the temperature. It is important to note, that A~(q) 
reflects the free energy density of different surface configura- 
tions of /~ but it is not a surface tension because no Goldstone 
modes, i.e., translational shifts without cost in free energy, ex- 
ist for /~ as they do for / + . Moreover, for q > one expects 
< Ao (q) ~ detE(g) for stability reasons. 



tension "f~{q) and an additional contribution that does 
not vanish for q — > and which depends parametrically 
only on the interaction potential w±2 between the two 
species and the planar density profiles d u p Ci (see Eq. (|p^|) 
for uj): 



A o {q) = -Q J {0,d,p Cl dnp C2 )+q 2 1 (q). 



(77) 



A~(q) describes the free energy required to deform the 
relative surface into a corrugated one with a wave- vector 
q in the presence of the microscopic interactions of the 
particles. For q = one has A~(0) > 0, which corre- 
sponds to the free energy needed to separate the flat equi- 
librium surfaces c\ and ci from each other. This is in ac- 
cordance with the facts, that A~(0) depends on w\i and 
significantly weakens for larger 8c\2 for which f% and /f 
decouple. In addition, at low temperatures A~(q) varies 
sensitively upon changes of 5c\2, but it hardly changes 
its character at higher temperatures. This can be ex- 
plained hcuristically by noting that for T < T c the dom- 
inant length scale is set by the diverging bulk correlation 
length £ so that the difference 8c\2 becomes irrelevant for 
the statistical weights. 
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A G (q) differs from 7 + (q) qualitatively (see Fig.[7|and 
note that A + (q,G = 0)/A+{q -> 0,G = 0) = 
7 + (g)/7 + (0) according to Ea. lf75|l ): for temperatures 
close to the triple point, it shows a monotonic increase 
implying that surface configurations / _ with nonzero 
wavelengths are energetically suppressed. But for higher 
temperatures a minimum at q =/= evolves. This mini- 
mum is also shifted towards longer wavelengths for fur- 
ther increased temperatures but does not change its 
depth. Thus, together with the behavior of 7 + (<j r ), this 
means that for low temperatures the mean surface / + is 
more easily excited thermally than the relative surface 
/~. But for higher temperatures, the thermal fluctua- 
tions have a stronger influence on /~ while / + becomes 
more rigid. This behavior is quantitatively controlled 
by the curvature corrections characterized by pn (see 
Eq. lj55l) ~) and thus Cn- The influence of Cn on A~(q) 
becomes mainly visible through a shift of the depth of 
the minimum, which increases strongly for larger values 
of Cm- Thus, A~(q) may even become negative for cer- 
tain values of Cn- In the absence of gravity, i.e. for 
G = 0, one has A~(q) — detE(g) and thus A~(q) < 
means that the system becomes unstable. By switch- 
ing off all interactions between the two species, i.e., for 
Wi2 = the simple expression 



A (9) = 9 



2 711(9)722(9) 
711(9) + 722(9) 



W'12 = . 



(78) 



emerges, so that from Eqs. 11741 and l|77|) the relation 
7~(9)7 + (9) = 711(9)722(9) follows. 

In the general case of nonzero W12 the expressions for 
7y in Eq. for 7l A 2 in Eq. (JHHJ, and for 7Y in Eq. (|B9|) 
together with 

Ti(9) := 7ii(9) + ^ + / 2(g) (79) 



and 



T 2 (9) := 722(9) + (80 ) 



lead to the following expression for 7 (q) as defined in 
Eq. GZJ (see Eq. ifglTjl for d, and Eq. (523 for #y) 



7 (9) 



cj(0,&p Cl d„p C2 ) -7(2(9) , T^q^q) 



7+(9) 



2d 1 d 2 J2W ij (q)+j- sym (q). 



The last term 7~ ym (q) (see Eq. <|C1|1 in Appendix[UJ) turns 
out to be the smallest contribution and it is determined 
by the contrast between the two species. The terms in 
Eq. 1)81(1 are listed according to their quantitative impor- 
tance. The behavior of 7~ (q) is determined mainly by the 



first and the second term, while the third one is about 
one order of magnitude smaller than the previous ones, 
and the last one may be smaller by even two orders of 
magnitude. 7^(9) captures the wavelength dependence 
of A~(q) (see Eq. {77|). A comparison between 7+ and 
7~ (see Fig-El) shows, that J~ (q) /"f~ (0) also exhibits a 
minimum at a nonzero wavevector but its depth increases 
with increasing temperature. Hence, for large values of 
Cjv, 7 _ (<z) an d even 7~(0), which depends on Cn, too 
(see Eq. (|81(lV may become negative which probably indi- 
cates a breakdown of the Gaussian approximation or even 
of the concept of a relative surface. Nevertheless, for in- 
creasing temperatures its minimum is shifted to smaller 
values of q, analogous to the behavior of 7 + (9)/7 + (0). 
Quantitatively, one finds 7^ (q) j 7" 1 " (q) < 0.1 for all values 
of q and temperatures T < 0.9 T c {l~{q)h + {q) < 0.2 for 
T < 0.99 T c ), so that r y~(q) is about one order of mag- 
nitude smaller than 7 + (?). Therefore, one may regard 
the mean surface / + to be more rigid than the relative 
surface / _ . 

Recently diffuse X-ray scattering data from the liquid- 
vapor interfaces of Bi:Ga, Tl:Ga, and Pb:Ga binary liq- 
uid alloys rich in Ga have been reported [13] . In order to 
interpret these data the authors put forward an expres- 
sion similar to Eq. lf?5|l (see Eqs. (8)-(ll) in Ref. |22|), in 
which, however, different than in Eqs. 1)75(1 and ((76(1 only 
the curvature correction profile pn 2 (f) of the segregated 
component was used without taking into account the pro- 
file ph 1 (u) of the majority component, i.e., ph ± {u) = 0. 
Thus it appears to be highly rewarding to reanalyze these 
experimental data in a future contribution on the basis 
of the present full statistical description. This descrip- 
tion might also provide an understanding of recent syn- 
chrotron X-ray reflectivity data on the interfacial width, 
broadened by capillary waves, of the liquid-liquid inter- 
face of nitrobenzene and water 0, 0|. However, this 
would require to extent the present analysis to binary 
dipolar fluids Furthermore, a recent analysis of the 
fluctuation spectrum of lipid bilayer shows similarities to 
our description in terms of two interfaces |3r| . These 
authors also define a mean and a relative surface (see 
Eqs. (5) and (6) in Ref. in order to take into ac- 
count conformations of the bilayer via fluctuation modes 
of the bilayer thickness and of the bending modes of the 
mean surface of the bilayer. Their choice of boundary 
conditions leads to a decoupling of these modes in real 
space and an effective free energy for the bilayer defor- 
mations on both short and long wavelengths is derived 
(Eq. (21) in Ref. |2fi|). The main difference to our ap- 
proach consists in their specific choice of boundary con- 
ditions for lipid bilayers, which cannot be applied to fluid 
interfaces. Accordingly, within the Gaussian approxima- 
tion, in our approach a mode decoupling between f + 
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Fig. 8: Normalized wavevector dependent surface free energy 
7~(g)/7~(0) as given by Eqs. 1771 and I81H for various tem- 
peratures T — 0.7 . . . 0.99 T c . The parameters are chosen as in 
Fig-El Similar to "f + (q), J~{q) shows a minimum which shifts 
towards smaller wavevectors but increases in depth upon in- 
creasing the temperature. Thus, the thermally activated cre- 
ation of additional surface area of the relative surface be- 
comes easier for longer wavelengths at higher temperatures. 
For large values of Cjv (see Eg. 1551 1 y~(a) becomes nega- 
tive which probably indicates a breakdown of the Gaussian 
approximation or even of the concept of the relative surface. 



In the limits G — > and q — > 0, one has f3Cij{q) = 
[GQ + + q 2 "f + (q)V^_ for all pairs i,j G {1,2} as already 
predicted in Ref. [37| (see also Eq. I|C10|) in Appendix[UJ| . 
However, our approach allows us to go beyond the limits 
G — ► and q — ► in order to obtain the interfacial struc- 
ture of a binary liquid mixture on smaller wavelengths. 
Although all expression derived above include the influ- 
ence of the external potential, i.e., gravity, we restrict 
our considerations in this section to G — in order to 
simplify the following discussion. 

One obtains from Eq. I|C7|) the positive function (see Eqs. 
for Tif ) 



1 eff A-(0) 1 [A o -(0)-g 2 7^( g )] 2 

/VCii(<?) 711 (q> q 2 q 2 K(0) + q 2 jg(q) 

(83) 

and a similar expression for 622(9) by interchanging 
the indices 1 «-> 2. As mentioned above, one has 
0q 2 Cu(qr^ -C 1) = 1/7+ (g). In the limiting case 
W12 = one obtains (iq 2 Cu(q) — l/ju(q). On the other 
hand, for Cu(q) (Eq. (JE2J) one has 



and /~ is achieved only in Fourier space, where / + is 
defined as a wavelength-dependent weighted sum of the 
Fourier modes of the surfaces fa and fa (see Eqs. <|71[l - 
(|73|l 1. Consequently, /+ consists of a sum of correspond- 
ing convolutions, which, in general, cannot be written as 
a sum of fa and fa with constant weights as it is done in 
Ref. . Nevertheless, in Subsec. lll Fl we have also intro- 
duced a concept similar to the one used in Ref. 36]. In 
Subsec. lTTTl the mean surface /* (see Fig .0 and Eq. |JS1}) 
is defined in real space in order to analyze Tih using the 
curvature expansion (Eq. I|43b|) 1 and in order to express 
TL in terms of fa and fa taking into account the full cou- 
pling between them. 

IV. CORRELATION FUNCTIONS 

From the diagonalization in Eq. I)73fl it is clear, that the 
surfaces / + and / _ are uncorrelated within the Gaussian 
approximation and thus statistically independent. In or- 
der to obtain insight into the structure of the original 
interfaces, one has to consider the correlation functions 



(/P(q)//(q'))=27r5(q + q')C ij ( g ), 



(82) 



where f3Cij(q) = E^ 1 (g) denote the matrix elements of 
the inverse matrix E _1 (q) (see Eq. I|59[) 1. The corre- 
sponding detailed expressions are given in AppendixICl 



[3q 2 C 12 (q,G = 0) 



Ao(Q)-g 2 7it(g) (g) 

7+(<z)A -(q) 



(84) 



however, Eq. 184fl does not exhibit the form of Eq. (|83|l . 
because 7*2 (q) , denned in Eq. (|C5|) , is also a positive 
function for all values of q. Thus, C\2{q) changes its sign 
at a certain value q , which depends crucially on Aj(0). 
This means that for q > q Q the Fourier modes of the 



surfaces fa and fa are anti-correlated (see Figs. 9 (c) and 



10 (c) I 



Figures 9 (a) - "Jcj show the correlation functions 



q Cij (q, G — 0) for different temperatures and the pa- 
rameter choices /ro = 1-002, Wo/w^ = 1.05, 
and Wo /wo = 0-5- Although the parameter differ- 
ences of the two components is small the correlation func- 
tions q 2 Cn(q) and q 2 C22(q) exhibit a different behavior 
for small wavelenths and temperatures close to the triple 
point (see Figs. 9 (a) and (b) I. This result indicates a 
structural difference between the surfaces ff and /| on 
short length scales which vanishes for higher tempera- 
tures. Similar the correlation function q 2 G\2{q) indicates 
a (anti-) correlation between the Fourier modes of the sur- 
faces for low temperatures which becomes weaker for high 
tempertures (see Fig. 9 (c) 1. 
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Fig. 9: Normalized height-height correlation functions 
q 2 Cij(q, G = 0) (Eqs. for various temperatures T = 

0. 7. . .0.99T C . The comparison between | (a)| and | (b)| reveals 
that the correlation of the Fourier modes of /f and f§ are sim- 
ilar at elevated temperatures but differ on short length scales, 

1. e., large q values, at temperatures close to the triple point 
due to differences between 711(5) and 722(17) (see Eg. 1761 1. 
For q > q (T) the correlation function q 2 Cu{q, G = 0) is neg- 
ative (see I (c)\ so that the Fourier modes of the two different 
surfaces are anti-correlated. Parameters are given in the main 
text. 



Fig. 10: Normalized height-height correlation functions 
q 2 dj(q, G = 0) of the Fourier modes / 8 c (q) for various concen- 
trations Xi of component 1 and for T = 0.7 T c . The choices 
for the interaction parameters are given in the main text. For 
high concentrations X\ > 0.6 one finds for q 2 Cu(q, G = 0) the 
well known structure for a height-height correlation function 
as predicted by classical capillary wave theory (see | (a)} . But 
for low concentrations X\ < 0.4 this function exhibits a peak 
for 2 < gr' 11 ' < 6 indicating the increasing influence of the 
second component on the interface /1 . This effect can be seen 
in reverse in | (b)| for small X2 = 1 — X±. Since Ci2(q) = C?2i(g) 
it does not depend strongly on the concentration (see | (c)| . 



Figures 10 (a) ■ TcT show the influence of the composi- 



16 



^3- 

■v. 



3- 



4.5 

4 
3.5 

3 
2.5 

2 
1.5 

1 

0.5 





1 

0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 



(a) 



X 1 = 0.6 






1 1 
Xj = 0.2 




Xj = 0.4 


— \ x '** 


Xj-0.6 " 




Xj - 0.8 _ 






(b) 


1 1 



(11) 

1 r o 



Fig. 11: Ratio of the correlation functions 6*22(9, G = 
0)/Cu{q,G = 0) and Ci 2 (q,G = 0)/C n (q,G = 0) for various 
concentrations X\ of component 1 and for T = 0.7 T c . Since 
all correlation functions attain the same value for q = 0, one 
has Cij(q = 0)/Cn(q = 0) = 1. One finds, that these ratios 
change significantly for different choices of the concentration: 
both, the slope for small q values and the limit q — > 00 are 
influenced by X\. Parameters are the same as in Fig. llUI 



tion X\ = 1 — X2 on the correlation functions q 2 Cij(q) 
at a fixed temperature T close to the triple point for the 
parameter choices To /r 1 ^ = 1.2, ajf 2 ' = 1.4, and 
Wo /wo — 0.586. We infer from Fig. 10 (b) that for 
low concentrations of species 1 the correlation function of 
the height of the interface associated with the more at- 
tractive component 2 does not show particular features at 
large q values whereas the interface of the less attractive 
component 1 seems to be more ordered by the stronger 
species. In Fig. 10 (a) this is indicated by a weaker decay 
of the corresponding height-height correlation function 
q 2 C\i{q) for X\ < 0.4. Since q 2 C\2{q) is a symmetric 
w.r.t. the label exchange 1 <-> 2, one would not expect 
a strong dependence of the corresponding height-height 
cross correlation function on the concentration. This ex- 
pectation is confirmed by Fig. 10 (c) 



V. SUMMARY 

We have considered the liquid-vapor or liquid-liquid 
interface region of a binary liquid mixture (Fig.^) which 
is characterized by two phase separating surfaces /i(R) 
for the two species i £ {1, 2} (Fig.^. We have obtained 
the following main results: 

(1) Based on a grand canonical density functional 
57[,oi(r), /02(f)] for binary liquid mixtures, we have de- 
fined an effective interface Hamiltonian 7i[/i(R), /2(R)] 
providing the statistical weight exp(— (3TL) for nonflat 
surface configurations (Eq. This approach takes 
into account and keeps track of both the presence of long- 
ranged dispersion forces (Eq. (|12|) '). smoothly varying in- 
trinsic density profiles (Fig.|3J), and the thermodynamic 
state of the system (Fig.^l. In particular, it captures 
liquid-vapor as well as liquid-liquid interfaces (see the re- 
marks (i)-(iii) in Sec.QJ. 

(2) Using a local normal coordinate system (Fig.01 
Eq. we have incorporated changes of the intrinsic 
density profiles caused by the curvatures of the fluctu- 
ating interfaces (Eq. (|43b|> U To this end we have intro- 
duced the concept of a mean surface /* (Subsec. lll Fl 
Fig.[3 and Eq. (|54fl ). Within this approximation the two 
surfaces fi and their width £j are used to form a nomi- 
nal surface /* for which the above mentioned coordinate 
change is applied without using additional parameters 
(see Sec.|H(iv)). 

(3) This approach leads to an explicit expression of 
"H[/i(R), /2(R)] in terms of the two surfaces (Ea. f2"K|) 
and Subsecs.HlEl and InF) . In particular it contains 
the coupling between the two surfaces based on the 
microscopic interactions between the two species (Ap- 
pendix lA"!) . 

(4) Within a Gaussian approximation the Hamilto- 
nian 7i[/i(R), / 2 (R)] takes a bilinear form (Eq. J5SJ)). 
In order to simplify the further discussion we define 
a mean surface / + and a relative surface /~ (Eqs. 
(|71|l and (JZH) which leads to a diagonalization of H 
(Eq. J73J)). Thus, the wavevector dependent free en- 
ergy density A + (q) (Eqs. (|74fl - H76|' ') of the mean surface 
/ + (Eq. 117111 1 generalizes the corresponding expression 
for the liquid-vapor interface of a one-component fluid. 
Therefore / + plays the role of the overall interface of the 
binary mixture which remains even in the special case 
fi = f2 = f + or f~ = 0, respectively, which is equiv- 
alent to a two-component system modelled by a single 
interface. A. + (q) contains a gravity part, GQ + (q), and a 
contribution stemming from the interactions among and 
between the species which can be considered as a wave- 
length dependent surface tension 7 + (q) (Eq. I|75|l . Ap- 
pendix^. 7 + (g) decreases as function of q, attains a 
minimum and increases again (Fig.[SJ). The minimum oc- 
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curs at smaller values of q and becomes less deep upon 
raising the temperature. This resembles the behavior of 
the wavelength dependent surface tension of the corre- 
sponding one-component fluid. The second g-dependent 
free energy density A~ (q) is linked to the relative sur- 
face /~. Even in the absence of the gravity (G = 0) 
and thus different from A + (q) it consists of an effective 
surface tension component "f~(q), and an additional con- 
stant contribution depending on the flat intrinsic density 
profiles and the interaction potential w\2 between the two 
species only (Eq. $771 ). This constant describes the cost 
in free energy for a separation of the two surfaces against 
the attraction between the two species. For temperatures 
close to the triple point A~(g) = A~(g, G = 0) increases 
monotonicly as a function of q, but for higher tempera- 
tures it developes a minimum that is gradually shifted to 
larger wavelengths (Fig.[7J|. 

The surface tension 7 _ (q) has a similar structure as 
7 + (g), but it develops a minimum the depth of which 
increases upon increasing temperature. Thus, depend- 
ing on the strength Cat of the influence of curvatures on 
the intrinsic profiles (Eq. JSEJO; 7~(<z) ma y even become 
negative signalling probably the breakdown of the Gaus- 
sian approximation or even of the mean interface concept 
for large values of Cm and certain temperatures (Fig.|HJ . 
However, the total energy density A~ (q) remains positive 
for all values of q. 

(5) Finally, we have discussed the Fourier transforms 
of the height-height correlation functions (Eqs. I|C7|I and 
ijnSJl). Figures 9 (a) - (c)| illust rate their temperature de- 
pendence whereas Figs. |10 (aj - Jc) demonstrate the influ- 
ence of the concentration on the height-height correlation 
functions (see Sec.H](v)). 



APPENDIX A: EXPLICIT FORM OF THE EF- 
FECTIVE INTERFACE HAMILTONIAN 

In this appendix we drop the tilde of Spf a (see 
Eq. Ij43b|l ^ and write 5p a instead and we frequently omit 
the full list of arguments an expression depends on. 
Hence, one should keep in mind that Sp a (S, u) and p c (u) 
depend on the set of normal coordinates. Here, we 
present our results for the effective interface Hamiltonian 
up to second order in the height displacements. Higher 
order terms and further details concerning the derivation 
can be found in Rcf. 31 1. 



pressed in terms of Agj = nijApj as 

2 

/ d 2 RH v (f(R)) = f / d 2 S J^Agj [Gf + Gt + Gf] . 

J A 1 J A j=1 

(Al) 

Similar as in Ref. [3j in order to proceed and for later 
purposes we define m 
.g 4 (S) = l + (V/ 4 (S)) : 
sions Ai = Ai(u, . . .) 



purposes we define moments with n > and the metric 
gi(S) = 1 + (V/i(S)) 2 (see Eq. gDJO of arbitrary expres- 



1 



and 



S n [Ai] ■= \fgj5n\Ai]- 



duu n A l (u,...) (A2) 



(A3) 



Without carrying out the curvature expansion we find up 
to second order in fy. 

Gf = - 2/? Si [d u p C] } + 2H j S 3 [d u p Cj } ( A4) 

- K j9j s4d aPci ]+ (/;) 2 

+ ({VfC) 2 + 2fC(2H 3 )^g-)6 2 [duP C] ], 



Q s = M + f c ZJoM 3 {2Hj) - §2[§pj] (A5) 



/9j V9j 
^HMSpA+AKj^UtPj] 



and 



- 1 



PW93 J-R 



:) 2 V/ J (S)V^(S, U ). 

(A6) 



Using the curvature expansion in Eq. 143 bfl for 5p up to 
quadratic order one obtains 



G*= 



/ 2Si\p H 2 
{2H, J f[—^~U 2 [PH ] 



<9i 



+ 



(/|(2^) 5 [ PH] ] + Kj 5i[pk,}) + 0{f 



(A7) 



By carrying out integration by parts with respect to the 
lateral coordinates one obtains for G^ in Eq. (|A1() 

l3 h[p Hl ]+0{f). (A8) 



G 



/9j 



where Afj means the Laplacian of the surface /j(S). 



1. Gravity Part 



Interaction Part 



After a transformation into appropriate normal coor- 
dinates the general form of the gravity parts can be ex- 



In this subsection we use the following notation: ex- 
pressions with an index i depend on u' and R' whereas 
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terms with an index j depend on u" and R". More- 
over we use the following symbolic notation: d' = d u >, 
V' = (ds^ , d s ' y ) , similarly d», V", and 

dS[p fi ,p fj } := d>p Ci d<>5pf. + d<5p fi d«p Cj + d<5p fi d«5p fj . 

(A9) 

In addition. we introduce the short notation 

^IMS'IMS")] = 4f[fi,fj] for * e {0,1,2} 

with (see Eq. for 7} and Eqs. (JT3J and for 

<■ ) 



,( fc )[ 



\fiJi] := wff(\T fi (S',u')-T f .(S",u")\) (A10) 



Wij[Scij] 

where wf^ = my. Similar to Eq. I|A2J) . it is convenient to 
use the following abbreviation for an arbitrary expression 
•Aij — -Aij (u , u ,...). 



w%\\S'-S"\,\u'-u" + 5c ij \), (All) 



Rm 
,(fc) 



dw'dw" w- fe) [/i, /,•] Aj {u',u",...) 



and similarly for [y^ij] for using w^- [<5cy] (Eq. l|All|l ) 
instead of tuff [f t , fj] in Eq. flAT0|) . uf ] still depends 

on S' and S"; hence, terms like Vs'oJf = V/oA (and 

similarly for V") are defined. The symbol (i/ ) <-» (j," ) 

< — > 

is introduced to shorten the formulae below. It states 
that the first part of the formula has to be repeated with 
interchanged labels, i.e., in each expression i is replaced 
by j and ' is replaced by ", and vice versa. 

After the transformation to normal coordinates the 
total expression for the interaction parts can be written 
as (Eq. J21) 



JJ d 2 R'd 2 R" H™(f(R,R')) 



■E 



d 2 S'd 2 S" 



w: 



OR 



(A13) 



with 



W- 



w 



- {uf[d.p^d,,p Ci \-w?>[d,p Ci d..p 0i ]) (A14) 

- 2HiU)f [u'd.p Ci d»p Cj ]+ K i wf[(u'fd.p Ci d..p Cj ] 
+ 2H i H j uf [u'u"d,p Ct d„p C] ] + ft/) «-» (j,") , 



: -o;f [dS[p fi ,p fj ]] (A15) 

2ff, W f [«' d8[p fi , p/,]] + Jfc 4 2) [(z/) 2 3<5[p /s , p/J] 
+ 2H l H 3 ojf [u'u" dS[p fi , p fj ]} + , 



and 



where 

w 0K 

y 



W 



(A16) 



,( 2 )r 



, 7 '[ftpcVz/fypj] - 2H i wf ) [u'd,p et V»6p J ] 



+ K l J 2) [{u l ) 2 d,p Ci V„5p. j \ 



(A17) 



Here, we already have omitted higher order terms which 
can be found in Ref. [3l| . Applying the curvature ex- 
pansion (Eq. (|43b|n for each density and its surface we 
obtain 



W 



as 



2Ht 



i ujf ] [d» p C] d> p Hi ] + Ki ojf [d» p Cj d> p K 



(2)r 



{2H l ) 2 [J 2 \d ll p C] d,p H ,] 
9", 

-Ki/)<-(j>") + 0(/ 3 ) 



u d>p Hi d»p Cj ] 



AHjHi ^^f ) [d / p Ht d"PH ] ] - uf\u' d> p Ci dn p Hj ] 



(A18) 



(A12) 

and up to second order 



w: 



an 



< > 

(A19) 



3. Hard Sphere Part 

In this subsection we use /* as introduced in Eq. 
and instead of h(p\, P 2) we write (see the introductory 
remarks in Subsec. fll Fl and Eq. 145(1 ) 



h{ Pl ,p 2 ) = h*{p*) 



(A20) 



and similarly hj = h*{fi),) and h* = h*(p*,). For the 
reasoning below it is not necessary to specify h*(p*) ex- 
plicitly. Using u = u* (distance from /*), S = S*, and 
(5p* = Sp*f* we find, after an integration by parts with 
respect to z, 



d 2 RH h (f(R,)) 



'A 

so that 



j d 2 R j dz (c* - z) [dh* f d zP } - dh* c d z p*] (A21a) 

J A J — OG 



d 2 RH h ^(R)) = / d 2 S Ap 



H a + FT 5 + R n 



(A21b) 

Similar to Eq. l(A2p we define the abbreviations (using 
Ap* and g* instead of Ap^ and <?.;) 



^l fc) [^*] 



(A22) 
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and apply the expansion 



dh} = dh*{ P * c + 5p*) » dh* + d 2 h* c 5p* 



d 3 K 



[Sp* 



(A23) 



to Eq. l|X2Tajl . With f® := f* - c* we find 

h 9 = -frf'M + a^W^M (A24) 

+ J® 2H* <pW [d uP * c ] - K*g* <p™ [dupl] 



1- 



^ ] [d u p* c 



FT 
and 



2^5 



, (A25) 



_ , duif® - 
Ap*# ./-ft*. V V5 



9/1} Vf*VSp* . 



(A26) 

Now we use the curvature expansion from Eq. I|43b|l for 
<5p*. Up to second order in /* one has 

H 4 = JL (2H*) [p^] + -L= X* rf) [ P y (A27) 



2 VF 



(2iT) 2 (- 



The contribution H n is treated similarly as (see 
Eq. I|A6[) ). The curvature expansion finally leads to 



H 



■R 



Af* (2H*)^ 1 1) [p* H } + 0((n 3 ).(A28) 



APPENDIX B: SECOND ORDER APPROXIMA- 
TION 

In this appendix we provide the explicit expressions 
up to second order which result from applying the equi- 
librium condition in Eq. (JJJ to the expressions given in 
Appendix^ Furthermore, the approximation 2H = 

9 " 2 (/,,; 1 + /J* + /„„(! + ID - -'./:,/,/,,) « A/ and 
the Gaufi-Bonnet theorem, J K(S)d 2 S = 2itxe, where 
Xe denotes the Euler characteristic of the surface /, are 
used. As a consequence, since we consider laterally flat, 
connected surfaces, we have xe — 0. Moreover it is 
convenient to present the terms in Fourier space using 
Eq. ll-'il'l) . Additional details are provided in Ref. • 



1. Gravity Part 

Together with the interaction terms we obtain from 
Eq. JXTJ with Eqs. IpOjl . (|A7|) . and (JXgJ by using the 



equilibrium condition in Eq. © several times 
G 



d 2 SH$({(S)) 



-in 



d 2 q 5>(<z)|/;(q)| 2 , (Bl) 
J'=l 



where 



Gj(q) := rrijApj ■ (l - 2q 2 5 [p H] ] 



(B2) 



Here, some expressions arising from the special treatment 
of the hard sphere part are not taken into account; we 
shall add them in Subsec. lB~3*l where the derivation of 
those terms is explained. The expressions in Eqs. I|B1|) 
and (|B2|) are identical to those derived previously for a 
single interface 



Interaction Part 



iv-j ! [q, 5cij] := w^\q,u' — u" + Scij) denotes the 
Fourier transformed interaction potential (see Eq. J23J), 
or the Fourier transformed integrals of Wij (see Eqs. <|13f) 
and PUP , respectively. Thus, similar to Eq. (|A12(1 for 
Aij ee Aij (u', u", . . .) we introduce 



du'du" [q, dc^] Aij{u', u", . . .) 

(B3) 



and 



S^ k \q,---)~^ k) (9,---)-^ k) (0,---), (B4) 

but we suppress the index c, because all quantities aA 
have been expanded in terms of / c and V/, respectively. 
Moreover, we have omitted the square brackets indicat- 
ing the functional dependence of u>M(q, [Aij(u', «",-..)]) 
(see also Eq. and we simplify . .) ee u;(. . .) 

due to w^ 0) (. ..) ee Wij (...). 

From Ea.fA13|l with Eas. l|AT4|) . (|AT8)l . and (|A"T9|) 

and by using the equilibrium condition in Eq. © we 
arrive at 

i2 q! i2 rill uG 



d'S'd'S" Hft.(f(S',S")) 



Wi 2 (g) V ai (g) 
:(<?) W 2 i(g) 



f (q) .(B5) 



The matrix elements and stem from different 
ui^(q, . . .), which include the planar density profile p c {u) 
and the first term ph(u) in the curvature expansion. In 
order to obtain a transparent presentation we define 



:=7«(g) + 7«(g) + 7«(g) 



and 



v«(«) 



7* 



(B6) 
(B7) 
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where (with the convention that quantities with an index 
i depend on u , while those with an index j depend on 
u"; see Subsec. lA~2)) 

7*0 (<?):= u(q,d<p Cz d„p Cj ) 
u^q^H.PHj) + 

7 X.(q) := -w(0,ap Ci a»p c ,) - 2q 2 u ! (0,p Hi d,p CJ ) 
+ 2q i GjW(0,p Hi d„p Cj ). 



(B8) 



and 



(B9) 



The constants Ky (Eqs. (|B6(I and l|B7|l ) are given ex- 



plicitly in Eq. (|B14(I and are introduced here for conve- 
nience in order simplify the calculations using = 7A, 

7! v 2 ^ 7 2 V i, 7* + 7^ = 9 2 7«, and + Ttf) = 9 2 ^ + 

(see Eqs. ifT^ll and JZHjl)- Since in 7^(9) the terms q 4 Kij 



are added with opposite sign, Eqs. I|B5|) - (|B7|) do not de- 
pend on . 

3. Hard Sphere Part 



Starting from Eq. (|A21b|) , we obtain up to terms sec- 
ond order in /* (except those which vanish identically 
due to the above-mentioned Gaufi-Bonnet theorem) 



d 2 RH h (f{H)) 



A 

Ap* 



d 2 S 



{f-c)^[d u pl 



(BIO) 



+ 2H* (^[d u p*} + cp^[p* H ]) - - (vrf^M 
+ (A/*) 2 (<pW\pU - <pU\p* B ] + \ <P?>[(P%?] 



Our aim is to express the right hand side of Eq. (|B10|) 
in terms of /j — Cj, i € {1,2}. To this end, we use an 
expansion of Eq. (|A21b|) around the equilibrium densities 
p Ci on the left hand side and an expansion around p* on 
the right hand side. A comparison of the terms leads to 
the relations 

2 

Y,djh(p Cl ,p C2 )Sp 3 = dh* c 5p* (Bll) 



£ d 2 3 h(p Cl , p C2 ) 5 Pl 8p 3 = d 2 h* c (8p*f . (B12) 



and 



Thus using the curvature expansion in Eq. I)43bjl for each 
Spi in Eq. I|B11|I one obtains similar equations which, 
for instance, relate the curvatures H 3 and H*. Using 



2Hi « A/,, 2H* » A/*, and /* - c* = (£ 2 /f 4 
£1 fi)/(£i + £2) from Eq. I|54|) one can express the con 



tributions in Eq. ijBlOfl including ip\ [p* H ] and Lp \p* H 2\ 
in terms of Af 3 . The same line of argument holds for 
Eq. HB12|I which yields 



]T K l3 A/, Af 3 » Ap* ^ [(pj,) 2 ] (A/*) 2 (B13) 
with the constants 

d 2 3 h(p Cl (u), p C2 (u)) pHt (u)ph s (u) , 

(B14) 



which were already used in Eqs. (|B8|) and (|B9|I and form 
the matrix K := (Ky 0i,je{i,2} ■ If we take into account 
the additional terms arising from using the equilibrium 
condition leading to Eqs. IjBll) and (|B5(I . we obtain for 
W£(f)=H£(f(S)): 



d 2 S?^(f) = ^- /d 2 9 ft(q) 



2q 2 



M{{q) 
Jft(q) 



2q 2 Bj2-^(q) + <i 4 - 

3=1 



f(q) (B15) 



with (see Eq. (JA3J)) 
J#?(g) := A Pl (6i [dihd u p Ci ] + q 2 S 1 [dihp Hi ]) , (B16) 



6 



and 



6+6 ' 

C?2 C?lC? 2 

did 2 d\ 



(B17) 



(B18) 



Furthermore, J4%(q) can be rewritten by using again 
the equilibrium condition Eq. |JSJ in order to express 
W^(f (S)) only in terms of the external potential V and 
the interaction potentials Wi 3 . Only the constant matrix 
K remains as a direct hard sphere contribution: 

2 

-J%{q)=G<g i {q)+Y,W i M (B19) 



where 



and 



:= miApj (<5 2 [<9 u p C! ] + q 2 S 2 [p Hi 



(B20) 



3%(g) ■.= u>M{0,u , d.p Ci d,. Pcj ) +q 2 ^ l \0 1 u'p Ht d, l p C] ) . 

(B21) 
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In combination with Eq. (|B2|) this gives the total gravity 
contribution in matrix form (see Eq. I|B18[1 ) 

G(g) := G ( Ql ( q) g^ q) )+^GBq 2 J2^(q) (B22) 



3=1 



with 



(B23) 



Along the same lines the interaction contributions can be 
expressed as (see Eqs. (IB6|) and flB7[ll 



W(g) 



with 



WM V 2 i{q) 
Via(g) W 2 i(<z) 



2g 2 © ^ #^-(g) (B24) 

*)3=1 



and similarly for E^ 1 (q) by interchanging the indices 1 
2, and 



Er 2 1 (9) 



-G 12 (g)+A-(0)- g 2 7l e 2 ff ( g ) 
detE 



(C8) 



Since all correlation functions share the same denomina- 
tor (Eqs. (IB2HI . GH), and (O), 



(C9) 



«r 1 Gn (9)7^(5) + 022(9)7^(9) 



detE(g) = detG(q)+A-(0)Gg+{q) 

'22 

-2Gi 2 (g)7if(?) 
+ g 2 7 + (<z)A -(g), 



one obtains for a vanishing external field strength, i.e., 
for G — > the following long-wave limit q — > for all 
pairs z, j e {1,2}: 



W«(?) := Wy (?) - 2g 2 ^ ^ tt (g) . (B25) 



fc=l 



APPENDIX C: EXPLICIT FORM OF THE COR- 
RELATION FUNCTIONS 

With Eqs. JEST), ©> and [jB"T7|l one has 



7alym(9) = ( Fife) - r 2 (g) - SW(q) ) (CI) 



with 



t+Gz) 

2 2 

<5^'(<?) == 2^ 2^ M (g)-2d2 ^#2,(<?). (C2) 
3=1 3=1 



In order to obtain compact expressions for the corre- 
lation functions (/i(q)/j(— q)) we use several abbrevia- 
tions. First, in order to simplify the matrix elements 
E-. 1 ^) which are related to the correlation functions (see 
Eq. JHU) ; we introduce the following expressions using 
Eqs.O-GSJ, I©-©) and (O (the argument q on 
the rhs is omitted): 



eff 



7iT(<7) := Tx-W-lS 

122 (?) ■= r 2 +sw- 7l e 



eff 



(C3) 
(C4) 



and 



^ : _ c-.-^)f^^) _ 7 -. (C5) 



This leads to the relations 



Etf(g) - G ij -(g) + (-l) i +U-(0) + ? 2 7^ f (g) (C6) 



resulting in 



, _ G 22 (g)+A-(0)+g 2 7 2 f(g) 
%1 (q) ~ detE 



(C7) 



G-fO 



1 



A+(g) 



+ 0(G,g 2 ), 



(CIO) 



which agrees with Ref. |37j . 

It is important to mention, that by using the sharp 
kink profile (Eq. (|26fl ) and by neglecting all curvature 
contributions at this stage, i.e., pn = 0, expressions for 
7 + and 7" arise which would not follow from the origi- 
nal sharp kink approximation introduced in Subsec TlI Dl 
Here, in this case one has #y = and Sfj = (see 
Eqs. 1C2() and jB2(ty ). which implies 7~ ym = and 
Gi k (q) — rriiApi. Furthermore, by using Swij(q, . . .) = 
Wij (g, . . .) — ibij (0, . . .) we obtain 



75 k (g) 



ApiApj 



8w l j(q,5c r 



r 



(Cll) 



and, with 7+< sk (g) = Ey 7// (9) 



, sk{q) = (Ap 1 ) 2 (Ap 2 ) 2 x 



g2 7 +,sk (?) 



(C12) 



Swn(q,0)8w 22 (q,0) - 8w\ 2 (q, Sc 12 



If all interaction potentials Wij have the same form and 

their amplitudes fulfill wi 12 ^ = \J roj 11 ' 'wf 2 ', all contri- 
butions to -f~< sk vanish with the only exception the non- 
vanishing contributions stemming from <5ci 2 and the dif- 
ference in the particle diameters. In addition, in this case 

with L(a) := *J (r^) 2 + ot 2 one finds a generalization of 
Eq. Jim or Eq. QJJ, respectively for Sc 12 ^ 0: 



i,3=l 



L 2 (Sc l3 ) 



( 



x 1 + 



g 2 L 2 ((5cy) r l0£( gL(Jc ii ) ) 



(C13) 
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